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Abstract. The critical behaviour of semi-infinite ci-dimensional systems with short-range 
interactions and an O(n) invariant Hamiltonian is investigated at an m-axial Lifshitz point 
with an isotropic wave-vector instability in an m-dimensional subspace of R d parallel to the 
surface. Continuum \<p\ 4 models representing the associated universality classes of surface 
critical behaviour are constructed. In the boundary parts of their Hamiltonians quadratic 
derivative terms (involving a dimensionless coupling constant A) must be included in addition 
to the familiar ones <x 4> 2 ■ Beyond one-loop order the infrared-stable fixed points describing 
the ordinary, special and extraordinary transitions ind = 4+ ^ — e dimensions (with e > 0) 
are located at A = A* = O(e). At second order in e, the surface critical exponents of both 
the ordinary and the special transitions start to deviate from their m = analogues. Results 
to order e' 2 are presented for the surface critical exponent /3J rd of the ordinary transition. The 
scaling dimension of the surface energy density is shown to be given exactly by d+m (0 — 1), 
where 8 = Vn/vi2 is the bulk anisotropy exponent. 
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Lifshitz points, i.e. multicritical points at which a disordered, a homogeneous ordered, 
and a modulated ordered phase meet, have been known since the end of the 1970s [1—4]. 
Appropriate n-vector </>| 4 -models representing universality classes of m-axial Lifshitz points 
were introduced at the same time; the simplest ones have a Hamiltonian Ti = J d d xCb(x) 
with density 

C b = | (A Q 0) 2 + I (V^) 2 + | (V Q 0) 2 + ^ + J |0| 4 , (1) 

where the position vector x = (x a , xp) has m- and (d— m) -dimensional components x a and 
xp, respectively, V Q and V/j denote the corresponding gradients and A a means the Laplacian 
V 2 

Although the possibility of studying the universality classes of these models in a 
systematic manner by means of expansions in d and m about general points on the line of 
upper critical dimensions d*(m) = 4 + m/2(0 < m < 8) had been realized already in 
1975 1 1 1, the enormous technical difficulties one encounters beyond one-loop order [5-8] had 
prevented a successful implementation of this programme until recently when a full two-loop 
renormalization group (RG) analysis was performed and the e = d* (m) — d expansions of all 
critical exponents were determined to order e 2 [9-13]. 
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Here we are concerned with the effects of surfaces on the critical behaviour at such m- 
axial Lifshitz points. The only previous studies of this problem we are aware of are restricted 
to the uniaxial (to = 1) Ising (n = 1) case and use either the mean-field approximation 
[13-15] or Monte Carlo simulations 1 16 1 for the ANNNI model. Let us consider semi-infinite 
systems with a boundary plane 03 at z — 0, where z > is the Cartesian coordinate along 
the inner normal on 03. Since x a and xp scale differently, two distinct basic orientations 
of the surface plane exist which we call parallel and perpendicular, depending on whether n 
is orthogonal to the a or the (3 subspace. We restrict ourselves here to the case of parallel 
surface orientation; the case of perpendicular orientation requires separate considerations and 
a distinct analysis 1171 . 

For semi-infinite |0| 4 models with an 0(n) symmetric Hamiltonian three distinct types 
of surface transitions occurring at the bulk critical point can be distinguished 1 18 19 1: the 
ordinary, special and extraordinary transitions. § Analogues of these surface transitions should 
exist also for the m-axial bulk Lifshitz points described by Hamiltonians with the bulk density 
(see footnote §). For the uniaxial Ising case TO = n= line? = 3 dimensions, Pleimling's 
Monte Carlo results 1 16 1 and the mean field analysis of 1151 lend support to this expectation. 
Our goal is to pave the ground for systematic field theory analyses of these transitions. 

To this end we need an appropriate semi-infinite extension of the bulk model with the 
density {0- For the short-range interaction case we are concerned with, it is justified to choose 
a Hamiltonian of the form (with W[ = M d_1 x [0, oo)) 

H= [ C b (x)dV + [ &(x)dA, (2) 

where C\{x) depends on 4>{x) and its derivatives. We must now (i) find out which 
contributions have to be retained in C\, (ii) determine the boundary conditions they imply, 
(iii) clarify the renormalization of the field theory and set up a RG approach in d* (m) — e 
dimensions, and (iv) derive the fixed-point structure, identifying potential fixed points 
describing the ordinary, special and extraordinary transitions. 

In the case of a critical point (corresponding to the choice to = 0), it is sufficient to 
include a term oc <f> 2 in C\\ other 0(n) (or I2) invariant contributions can be shown to be 
redundant or irrelevant 1 1 81 1191 . However, in the (m ^ 0) case of a Lifshitz point, this is not 
sufficient; we must take 

£ 1 (x) = ^cj> 2 + ^(V a <j>) 2 . (3) 

Power counting tells us that Act -1 / 2 is dimensionless. Hence it is scale invariant at the 
Gaussian fixed point and potentially infrared relevant for e > 0. All other contributions, 
notably terms oc (Vpcp) 2 , oc 4>W a <p or oc can be ruled out by symmetry or shown 

to be irrelevant or redundant 1 17]. The field theory defined by equations Q-Q satisfies the 
boundary conditions (valid in an operator sense l fl8lll9l ) 

d n <j>= (c-AA a )0. (4) 

This carries over to the free propagator G(x, x'), whose Fourier transform, G, with respect to 
the d — 1 coordinates parallel to the surface reads, in the disordered phase, 



G{p;z,z') = ^- 



(5) 



ip\z-z'\ _ C +MPa\ K P c -k v (z+z') 
C+\\p a \ 2 +K p 

§ The occurrence of the extraordinary and special transitions requires that the surface dimension d — 1 is sufficiently 
high so that long-range surface order is possible in the presence of a disordered bulk. 
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with 

k p = ^T + p\p a \2 + \ P/3 \ 2 + a\p a \ 4 , (6) 

where p a is the m-dimensional a component of the wave-vector p E K d_1 . The back 
transform of the part depending on \z — z'\ is the free bulk propagator G b (x — x'). At 
the Gaussian Lifshitz point f = p = ii = 0, it takes the scaling form 

G b (x) = \ X/3 \- 2+e a- m ^^ m ,ci(^ 1/4 \^\\x \- 1/2 ) . (7) 

Here the scaling function $ m ,d (v) is a generalization of a generalized hypergeometric 
function (a Fox- Wright i?/>i function, cf equations (10)— (13) of [10|). In the special cases 
c = A = 0orc^ooat arbitrary A > 0, G(x,x') reduces to the Neumann or Dirichlet 
propagator, respectively, whose z + z' dependent parts reduce to ±Gb (x — x' + 2z' n). 

Utilizing these results, and employing dimensional regularization in conjunction with 
minimal subtraction of poles, we have performed a two-loop RG analysis of the model (1)- 
(3) in d*(m) — e dimension. 

Its main results are as follows. To renormalize the multi-point correlation functions 
q(n,m) involving N fields </> off and M fields s = cj>{x e *B) on the boundary, the 'bulk' 
re-parameterizations known from I9I I10I . 

= 0ren ) <J = Z a (J , T - f LP = fl 2 Z T T , 

ift-pLp) °~ 1/2 =»Z p p, u&- m ^F m<E = pfZ v u, (8) 
must be complemented by 'surface' re-parameterizations of the form 
cjP = {Z^Zx) 1 ' 2 <j>f en , c-c sp = fiZ c c, Xa- 1 ' 2 = X + P x (u,X,e) , (9) 
where the surface renormalization factors Z\ and Z c depend on u and X. The function 

oo oo oo 

Px(u, X,e) = P\'~ J) W e " j =J2J2 P { x h ~ r ' k) u* X k (10) 

does not vanish at A = 0; although the one-loop coefficient Pj^' l ' (A) vanishes at A = 0, the 

graph o °f (<P 4> 3 ) yields a non-zero pj 2, 1,0 \ Thus a contribution oc (V Q 0) 2 

to C\ gets generated under the RG even if it was originally absent. 

The fixed points V* ld , P* p and P* x describing respectively the ordinary, special, and 
extraordinary transitions must lie in the cA plane at [r,p,u) — (0,0, u*), where u* is the 
nontrivial root of the bulk beta function f3 u (u) = fidfj,\ u, computed to order 0(e 2 ) in 1 10 1. 
For u — u*, the beta function (3\(u, A) = fid^\ Q X turns out to have an infrared-stable root at 

A* = -2ePf .-i:°)/p(i.-i!i) + (e 2 ) , (11) 

with a correction-to-scaling exponent 

oj x = (0 A /3 A )(«*, X*) = -pf" 1 ^ u * + 0(e 2 ) = e + 0(e 2 ) . (12) 

n + 8 

Thus P* rd , V* p , and P* x are the fixed points at c = oo, 0, and — oo displayed in figure^ 

Upon exploiting the RG equations implied by the above re-parameterizations and (|9j 
in a standard fashion, one concludes that the critical surface exponents of the special transition 
can be expressed in terms of bulk exponents and 

A[^} = (d-m-2 + r lL2 + 7 ? *< sp + m0)/2 = ff/uw (13) 
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A* = 0(e) 



1+ C 



Figure 1. Schematic picture of the RG flow in the cA plane if m > 0, showing the fixed points 



and 



A[s»] 



2 + m6- r)*- sp 



(14) 



the scaling exponents of (fy^ (x) and the boundary energy density e*(x) = [0 ffl ' (x)} 2 / % 
respectively, where the superscript '*,sp' means values taken at V* p . The e expansions of 
these exponents, like those of the bulk exponents, turn out to be independent of m to order 



e, but can be shown to be m dependent at 0(e ) 1171 . Thus, rj^ 



n+2 
' n+8 



e + 0(e 2 ) 



rr p + 0(e 2 ).|| 

These statements about the m-dependence apply equally well to the surface critical 
exponents of the ordinary transition. To demonstrate this via explicit 0(e 2 ) results, note 
that /?J' rd can be expressed quite generally in terms of standard bulk exponents v^, ?7l2 (or 
(3i), and a single additional anomalous dimension rj* ^ as 



0T d = (vL2/2)(d - m + n L2 + m6 + r)l <00 ) 
Our two-loop result for T)1 oo I s 
n + 2 



0L + Wl (1 + Vtoo /2) 



'/l. 
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{l + u* [i!(m) - J»(m)] } + 0[(u*) 



(15) 



(16) 



Here u* is the fixed-point value whose e expansion to 0(e 2 ) is given in equation (60) of 
1 10 1 while J u (rn), defined by equations (49) and (50) of that reference, is one of the four 
single integrals (j^, j a , j p , J u ) in terms of which the two-loop series coefficients of the bulk 
exponents were written there [see its equations (43)-(45) and (50)]. Finally, 

2 10+m 7r 6 + 3m / 4 r(™/2) 



dv i 



1 ; 



Jo 



(17) 



r(2-m/4)r 2 (m/4) Jo 

is a similar new integral which can be reduced to a single one. (Upon rewriting ^dv f Q dy as 
Jg 00 dy f^°dv, the latter v integration can be performed analytically to obtain a combination 
of hypergeometric functions.) 

2-! „„h ^ _ l _ n+2 i ^ (e 2 ) of /3f and the surface crossover 



|j The implied O(e) results = \ - f + 0(e 2 ) and <S> = § 
exponent <E> may be gleaned from equations (3.156e,b) of 1 18 1. 



n+2 _£ 
n+8 4 " 
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Combining these results yields 

n + 8 16(n + 8)' :i ^ Lm + 2 J 

+ 96(8 + n) \ji{m) - J„(m) Je 2 + 0(e 3 ) . (18) 

Just as for the above-metioned four integrals of 1101 . the values of ji(m) can be 
computed analytically for the special choices m = 2, m — 6 and m — > 0. This yields 

, , 1 . . In 3 , , 2 27 

ii(0) = -, ji(2) = 1 _ — , ^(6) = -- + 2 In- . (19) 

To determine j\ (m) for other values of m we had to resort to numerical means of the kind 
utilized in |10|. For the uniaxial case m = 1, we obtained ji(l) = 0.47289(1). Note also that 
in the limit m — > 0, the result Jl 8i reduces to the familiar one for the standard semi-infinite 
\(f>\ 4 model, given in equation (IV.35) of 1 20 1 . 

Let us briefly explain how the above results were obtained. Since the fixed point V* rd 
is located at c = oo, the ordinary transition can be investigated without having to retain the 
full dependence on c and A. To see this, note that the free propagator and the regularized 
bare G^ N ' M ^ become independent of A in the limit c — > oo, and satisfy a Dirichlet boundary 
condition, which carries over to the renormalized theory. The long-scale behaviour of the 
Q( N > M ) w ith a nonzero number of dr 5 can be inferred from the theory with c = oo and 
A = via the near-boundary behaviour of the operator <j>. To this end, one considers 
correlation functions involving an arbitrary number of the operators <p and d n <f>, and then 
uses the boundary operator expansion (BOE) 4> lcn (x<s + zri) « C or d(z) d n cf> lcn (x<s)- The 

renormalized theory requires in addition to the bulk re-parameterizations, the multiplicative 
re-parameterization 0, 

d n 4> = [z ll00 («) z^u)] 1 ' 2 a n ren (20) 

and an additive surface counter-term subtracting the primitive divergence (oc p 2 ) of 
(d n 4>d n <j>). The resulting RG equations imply scaling and yield the behaviour C or d(z) ~ 
z 1+Vl ,°=/ 2 , where 77*00 is the fixed-point value of the exponent function associated with 
Zi iOC . A straightforward consequence is that the exponents characterizing the leading infrared 
singularities of the G^'^ can be expressed in terms of (4 independent) bulk critical indices 
and a single surface one, namely, rjl ^ or (3° ld . Upon computing Zi )O0 and its exponent 
function to two-loop order and making extensive use of the results of |9| and 1 10], we arrived 
at equations J 1 5I > — d 1 81. 

Let us also note that the scaling dimension of the surface energy density at the ordinary 
fixed point is given exactly by 

A old [e <s }= d + m(6 -1) . (21) 

That is , the leading thermal singularity of e 23 has the bulk-free energy form ~ |t| 2 ~ Ql with 
Q!£ = Vh2(d — m + m9). The result can be obtained in a variety of ways, namely: (i) by 
generalizing the analysis given in appendix C of 11211 . (ii) by showing that the operator with 
smallest scaling dimension appearing in the BOE of the energy density e(x) is the component 
T zz of the stress-energy tensor [whose scaling dimension is given by equation (1211 1 and (iii) 
by proceeding as in the derivation for the m = case given in section III.B of |22 1. 

The results dl6> and Jl 8i can be combined with known bulk results [ 1 1 to estimate 
the values of surface critical exponents such as /3° rd for d = 3. In the uniaxial Ising case 
m = n = 1, equation Jl 81 becomes 

rjl ^(m = n = 1) = -0.3333e - 0.1804 e 2 + 0(e 3 ) , (22) 
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which gives rj'l ^ ~ —0.906 if we set e = 3/2 (i.e. d = 3), truncating the series at order 
e 2 . Inserting this value together with the bulk estimates (3^ ~ 0.246, — 0.746 and rj^ 2 — 
0.124 of 1 10 1 into equation dl5l and the analogous expressions 7° rd = (1— Tjh2 — Vi oo/2) 
and 7°J d = — ^L2 (1 + ?7L2 +i)i a! ) for the surface susceptibility exponents yields 

/3° rd ~ 0.65 , 7 f d ~ 0.99 , 7 ^ d ~ -0.2 (m = n = 1, d = 3) . (23) 

Owing to the low order e 2 of the available series expansions and the large value e = 3/2 
involved, these estimates cannot be trusted to be very precise. (They inherit, in particular, 
any uncertainty of the inserted bulk exponents.) However, they compare reasonably well 
with Pleimling's recent Monte Carlo estimates OH /3f d = 0.687(5), j? ld = 0.82(4) and 
7 ™ d = -0.29(6). 

In summary, we have identified the continuum models that represent the universality 
classes of the considered ordinary, special and extraordinary (surface) transitions at m-axial 
bulk Lifshitz points, clarified their fixed point structure and presented two-loop RG results. A 
more detailed account of this work will be presented elsewhere 1171 . 
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